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Abstract 

Several research groups are giving special attention to quantum walks recently, be- 
cause this research area have been used with success in the development of new 
efficient quantum algorithms. A general simulator of quantum walks is very impor- 
tant for the development of this area, since it allows the researchers to focus on the 
mathematical and physical aspects of the research instead of deviating the efforts 
to the implementation of specific numerical simulations. In this paper we present 
QWalk, a quantum walk simulator for one- and two-dimensional lattices. Finite 
two-dimensional lattices with generic topologies can be used. Decoherence can be 
simulated by performing measurements or by breaking links of the lattice. We use 
examples to explain the usage of the software and to show some recent results of 
the literature that are easily reproduced by the simulator. 
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C compiler. Successful tests were performed in Linux and Windows. 

RAM: Less than 10 MB were required for a two-dimensional lattice of size 201 x 201. 

About 400 MB, for a two-dimensional lattice of size 1601 x 1601. 

Supplementary material: Several examples of input files are provided. 

Keywords: Quantum walk; quantum computing; Quantum Mechanics; simulation; 

double-slit; broken links; C. 

PACS: 03.67.Lx, 05.40.Fb, 03.65. Yz 

Classification: 

Subprograms used: The simulator generates gnuplot scripts that can be used to make 
graphics of the output data. 

Nature of problem: Classical simulation of discrete quantum walks in one- and two- 
dimensional lattices. 

Solution method: Iterative approach without explicit representation of evolution 
operator. 

Restrictions: The available amount of RAM memory imposes a limit on the size 
of the simulations. 

Unusual features: The software provides an easy way of simulating decoherence 
through detectors or random broken links. In the two-dimensional simulations it 
also allows the definition of permanent broken links, besides calculation of total 
variation distance (from the uniform and from an approximate stationary distribu- 
tion) and the choice between two different physical lattices. It also provides an easy 
way of performing measurements on specific sites of the 2D lattice and the analysis 
of observation screens. In one-dimensional simulations it allows the choice between 
three different lattices. Both one- and two-dimensional simulations facilitates the 
generation of graphics by automatically generating gnuplot scrips. 

Additional comments: An earlier version of QWalk was first presented in [1]. 

Running time: The simulation of 100 steps for a two-dimensional lattice of size 
201 X 201 took less than 2 seconds on a Pentium IV 2.6GHz with 512MB of RAM 
memory, 512KB of cache memory and under Linux. It also took about 15 minutes 
for a lattice of size 1601 x 1601 on the same computer. Optimization option -02 
was used during compilation for these tests. 

References: 

[1] Marquezino, F.L. and Portugal, R., QWalk: Simulador de Caminhadas Quanticas, 
in Proceedings of 2nd WECIQ, pages 123-132, Campina Grande, Brazil, 2007, 
IQuanta. 
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LONG WRITE-UP 



1 Introduction 

A very successful approach used to solve some problems in classical computing 
consists in using algorithms based on random walks. In fact, many computa- 
tional problems are solved more efficiently by this technique [1]. In the 1990s 
the model of discrete quantum walks has been developed by Aharonov et al [2] 
and the continuous model has been developed by Farhi and Gutmann [3] . The 
good results obtained with random walks in classical computing motivate us 
to investigate the applications of quantum walks to the development of new 
efficient quantum algorithms. 

Quantum algorithms based on quantum walks have already been developed 
with great success. Some examples are the search algorithm by Shenvi et al [4] 
and the algorithm for element distinctness by Ambainis [5], both using the 
discrete-time model. Kempe [6] showed that a discrete-time quantum walker 
can cross a hypercube of dimension d>3 exponentially faster than a classical 
random walk. Farhi et al [7] have developed an efficient quantum algorithm 
for the Hamiltonian NAND tree based on the continuous-time quantum walk. 

Together with amplitude amplification and Fourier transform, quantum walks 
are one of the strategies for developing efficient quantum algorithms. Algo- 
rithms based on quantum walks are better than the Grover algorithm in some 
cases [5]. Apart from the applications on the development of efficient quan- 
tum algorithms, which are important from a Computer Science perspective, 
the quantum walks have also properties that justify their investigation from 
the point of view of Physics. 

A general simulator of quantum walks is important for the development of 
this research area. Without such a simulator, the effort of the researchers 
are deviated to the implementation of specific numerical simulations while it 
should the focused in the physical and mathematical aspects of the research. 
The QWalk simulator allows the scientific community to perform important 
simulations on quantum walks with simple commands and even facilitates the 
generation of plots to visualize the results. In its present form, the QWalk 
simulator can reproduce most of the simulations present in research papers. 

In Sec. 2 we briefly review the discrete quantum walks and some related defi- 
nitions. In Sec. 3 we present the QWalk simulator, using examples from recent 
results of literature in order to explain its usage. In Appendix A we describe 
some options of the simulator that were not used in the examples. 
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2 Quantum walks 



In the discrete one-dimensional quantum walk we consider a free particle — a 
walker — moving at each time step along a one-dimensional lattice. The direc- 
tion of movement is given by an additional degree of freedom, the chirality 
of the particle, which can take two values, analogous to the result of the coin 
tossing in the classical random walk. In the two-dimensional quantum walk, 
similarly, the walker moves at each time step along a two-dimensional lattice. 
In this case two additional degrees of freedom are necessary, so that one may 
decide among four kinds of movements. In this Section we describe in details 
only the two-dimensional case. The one-dimensional walk, nevertheless, can 
be seen as a simple particularization of the case here exposed. 

The Hilbert space considered in the two-dimensional walk is 7^2 ® "^2 ® "^oo, 
where 7^2 ® 7^2 is the Hilbert subspace associated with the chirality of the 
walker, i.e., the coin, and Tioo is the Hilbert space associated with the position 
of the particle over the lattice. The basis for the coin subspace is Be = {|j, k) : 
J, k e {0, 1}} and the basis for the position subspace is Bs = {|m, n) : m, n e 
Z}. 

The generic state of the quantum walker at time t is 

1 oo 

l*W)=E E ^3,k;m,n{t)\j.k)\m,n), (1) 

j,k=0 m,n=-co 

with i]jj^k-rn,n{t) G C and Ei,fc Em,n l'0j,fe;m,n(^)P = 1- The cvolutiou of the 
system over time is given by a unitary operator U = S o {C ® Ip), where S 
is the shift operator, Ip is the identity operator which acts on the position 
subspace and C is the coin operator, which acts on the 7^2 Ti.2 subspace. 

In this paper we address two different shift operators for the two-dimensional 
walk. The first of them has already been described in [8] and is given by 

1 +00 

^a-E E \j,k){j,k\^\m+i-iy,n + i-l)'){m,n\. (2) 

j,k=Q m,n=-oo 

One may notice that according to this shift operator the walker always moves 
along the diagonals of the mathematical lattice — it moves along the main 
diagonal if the value of the coin is either |00) or |11), and along the secondary 
diagonal if the value of the coin is either |01) or |10). In this case we say that 
the physical lattice is diagonal (in relation to the mathematical lattice). 

We may also consider the possibility that, at time t, a link from site (m, n) to 
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(a) Diagonal lattice 



(b) Natural lattice 



Fig. 1. Part of the lattice for a two-dimensional quantum walk, showing a broken 
link. 

one of its neighbours is broken. Then we need the functions 

^ , , , I (— l)-' if hnk to site m+ (—l)'',n + (—1)'^ is closed, , , 
Ci{j, k; 171,71) — < (3) 
I if link is open. 



and 
^(j, A;;m,n) 







-1)*^ if link to site m + (-!)■', n + {-!)'' is closed, 
if link is open. 



(4) 



with j, k e {0, 1}. Whenever - j,l - k;m + (-1)-^, n + (-1)'') = we 
must impose Ci{j, k; m, n) = 0, and similarly for £2- The technique of broken 
links for quantum walks was developed originally by Romanelli et al [9] and 
later generalized for the two-dimensional case by Oliveira et al [8]. 

In Fig. 1(a) we have part of the lattice used in the two-dimensional quantum 
walk, with shift operator Sa- The mathematical lattice is depicted by a dashed 
line and the physical lattice is depicted by a full line. In the example we have 
a broken link between sites (m, n) and (m -|- 1, n -|- 1). 

If we apply the evolution operator Ua to state (1) and include functions Ci 
and £2 as in [8] , we obtain the evolution equation 

i'l-j,l-k;m,n{t + 1) = 



j',fe'=0 



We can also define a second shift operator in order to obtain a physical lattice 
that coincides with the mathematical lattice. This operator is given by 



Sb=Y. E \3A){hd\® m + {-iy{l-5j,d),n + {-iy5j,d){m,n\. 

j,d=0 m,n=-oo 



(6) 
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Wc will observe through the examples of Section 3 that the probability dis- 
tributions obtained with operator Sb differ from those obtained with operator 
Sa only by a rotation of 7r/4. It will also be clear later on that, altough the 
movement described by Sb may be more intuitive, operator Sa has also some 
advantages. 

If we intend to include the possibility of broken links in this second lattice, we 
need the function 



d; m, n) 



-ly if link to m + (— 1)-'(1 — ^j^d), n + {—lySj^d is closed, 
if link is open, 

(7) 

with j,d e {0,1}. As in the previous case, we must impose C{j,d]m,n) = 
whenever C{1 - j, 1 - d; m + (-1)^(1 - 5,-,^), n + = 0. 

If we apply the evolution operator Ub to state (1) and include function £, we 
obtain the evolution equation 

'4'l-j,l-d;m,n{t + 1) = 
1 

j',d'=0 

where ® is addition modulo 2. 

In Fig. 1(b) we have part of the lattice used in the two-dimensional quantum 
walk with shift operator Sb- In this example we have a broken link between 
sites (m, n) and (m, n -|- 1). 

For the one-dimensional simulations, Qwalk provides three kinds of shift op- 
erators, which defines three different lattices: the infinite line, the finite one- 
dimensional lattice with reflecting boundaries and the cycle [10,11]. 



2.1 Average distribution and mixing time 



By suitably defining permanent broken links we may study the quantum walk 
on a finite lattice with interesting topologies, such as a square box with re- 
flecting boundaries [12]. Let 

P{m,n,T)= E \i^j,k;m,n{T)\' (9) 
j,k=0 

be the probability of finding the walker on the site (m, n) of the lattice at 
time T. In a classical random walk on a box this probability distribution 
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converges to a stationary distribution. On the other hand, since the quan- 
tum walk must preserve unitarity, it does not present this convergence prop- 
erty. However, defining the average probabihty distribution as P(m, n, T) = 
^ 121=0 P{fn,n,t), Aharonov et al [13] have proved that P{m,n,T) converges 
in the hmit T ^ oo. Hence, define 



7r{m,n) = hm P{m,n,T). (10) 

T — ^oo 

The rate at which the probabihty distribution approaches the hmiting distri- 
bution 7r(x) is captured by 

Definition 1 (Mixing time) The mixing time of a quantum Markov 
chain is 

M, = min{T|Vt > T, \\Pt - n\\ < e}, 
where \\A—B\\ = J2mn \A{m,n)—B{m,n)\ denotes the total variation distance. 

Usually depends on the initial conditions. 



3 Using the software 



The QWalk simulator is quite easy to install on a Linux-like environment. 
One just needs to download the source codep] uncompress it in any directory 
and finally use the command make. The source code of the simulator was 
also successfully compiled in the Microsoft Windows operating system, using 
Dev-C-|--|- 4.9.9.2 compiler, which can be downloaded for free in Internelp] 
Detailed information on how to compile the simulator, as well as pre-compiled 
versions of it, are also provided in the web-site. Together with the downloaded 
files the user with knowledge in C programming also finds information on how 
to change the source code. 

The simulator consists of three tools: qwld simulates quantum walks in one- 
dimensional lattices; qw2d, in two-dimensional lattices; and qwamplify im- 
proves the visualization of the plots generated by qw2d by amplifying some 
regions. In this Section we give some examples from recent results of literature 
in order to show how to use the simulator. Other examples are also provided 
with the downloaded files. 



^ http://qubit.lncc.br/qwalk 
^ http://www.bloodshed.net. 
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3.1 Simulations in two-dimensional lattices 



In order to use qw2d it is necessary to write an input file in any ASCII text 
editor. This input file consists of keywords that define the simulation options. 
Most important keywords are explained in the examples of this Section, and 
the ones not covered here are discussed in Appendix A. 

After creating the input file, say file, in, one just needs to type qw2d f ile . in. 
The results are stored in some output files. The f ile .dat output file contains 
the final probability distribution. The file-wave . dat output file contains the 
final complex amplitudes of the wave function. The f ile-pb.dat output file 
contains the approximate stationary distribution, when it is requested. The 
file-screen.dat output file contains the data observed in the screen, when 
it is requested. The file.sta output file contains certain statistics such as 
variance, standard deviation, average and total variation distance from an 
approximate stationary distribution and from the uniform distribution. The 
file. pit output file is a gnuplot script. The postscript files it generates de- 
pend on the options used. Typically these files are: f ile-3d. eps, the 3D plot; 
f ile-2d . eps, the contour plot; file-screen . eps, the observation screen pat- 
tern; f ile-pb . eps, the approximate stationary distribution plot. Further plots 
can also be manually generated by the user with some knowledge on gnuplot 
or any other similar tool. 



3.2 Simulation of a double-slit experiment 

In Fig. 2 we have the result of a simulation of the double-slit experiment 
with quantum walks, reproducing some results recently obtained by Oliveira 
et al [14]. This simulation took less than 2s on a Pentium IV 2.6GHz with 
512MB of RAM memory, 512KB of cache memory and under Linux. In order 
to perform this simulation with qw2d the input file must have the following 
keywords, all in uppercase: 



It is not necessary to place the commands exactly as in the example. The 
user may skip lines between different commands or even write everything in 
a single line. It is important, however, to keep all the keywords in the main 
section of the input file, delimited by a BEGIN keyword and an END keyword. 
Otherwise, the keywords will be interpreted as comments. The plots of Fig. 2 



BEGIN 
COIN HADAMARD 
STATE HADAMARD 
STEPS 100 

END 



BLPERMANENT 

SCREEN 60 -100 60 100 

LATTYPE DIAGONAL 



8 



100 




X 



(a) 3D plot (b) Contour plot 

Fig. 2. Probability distribution after a double-slit experiment. An amplification 
factor 5 was used for x > 20, in order to improve visualization. 

were generated by gnuplot, with a script provided as output of qw2d. The 
generation of the plots took about 30s. 

The COIN keyword defines the coin used in the simulation. In the previous 
example we chose the Hadamard coin. We could also have chosen the Fourier 
or Grover coins, with the options FOURIER and GROVER, or even an arbitrary 
unitary coin. In the latter case, the option to be used is CUSTOM. An extra 
section in the input file in also required in this case in order to define the 
arbitrary coin, see Appendix A. 

Analogously, the STATE keyword defines the initial state of the simulation. In 
the previous example we chose the state that provides a better spreading for 
the Hadamard coin. We could also have chosen the corresponding states for 
Grover or Fourier coins, or even an arbitrary initial state. 

The STEPS keyword defines the number of iterations that will be carried out 
by the simulation. In the previous example the walker performed a hundred 
simulation steps. The user should recall that the larger is the simulated time 
without boundaries the farther from the origin the particle can be found, 
which means that a bigger space in computer's memory must be reserved for 
the simulation. Therefore, this keyword may increase not only the running 
time but also the memory requirements of the simulation. It will be explained 
later on how to avoid this memory consume when we fix certain boundaries 
in the lattice. 

The slits were created with help of BLPERMANENT keyword in the main section of 
the input file. With this command we declare that some links in the simulation 
will be permanently broken. In order to define the position of those broken 
links we use a separated section in the same input file, with the commands 
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BEGINBL 

LINE 20 100 20 7 

LINE 20 5 20 -5 

LINE 20 -7 20 -100 
ENDBL 

The command LINE xO yO xl yl isolates all the points that pass through the 
segment from {xO,yO) to {xl,yl). The line of isolated points may be parallel 
to the axes X or y, or may have an angle of 45 degrees with one of those. It is 
also possible to isolate a single point with the command POINT xO yO. With 
combinations of these two commands one can simulate not only single- and 
double-slit experiments but also the evolution of the walker with an enormous 
variety of boundaries. In this first example we have a slit at (20, 6) and another 
at (20,-6). 

An observation screen may be defined with the SCREEN keyword, which must 
be followed by the initial and final coordinates. The screen may be defined 
parallel to axes x or y, or with an angle of 45 degrees with one of those. In 
the example the screen goes from (60, —100) to (60, 100). 

Since the fraction of the wave that passed through the shts in the previous 
example was very small, the visualization obtained in the first place by the 
simulation turned out to be quite precarious. In order to solve that, we used the 
qwamplify tool to amplify by a factor of 5 the whole region where ,t > 20. This 
tool may be used by typing something like qwamplify file . dat [options] . 
The software creates a safety copy of file. dat and replaces it by a new one, 
with the appropriate part of the wave function amplified. Help concerning the 
available options can be obtained simply by typing qwamplify. 

The LATTYPE DIAGONAL command declares that the shift operator Sa is to be 
used, see Eq. (2). The shift operator Sf,, from Eq. (6), is the default but can 
also be exphcitly declared with the LATTYPE NATURAL command. The natural 
lattice, when compared to the diagonal lattice, gives probability distributions 
rotated by an angle of 45 degrees. We can note this behaviour in Fig. 3, 
the simulation of the Hadamard walk for a hundred steps without slits. Note 
that there is a region on the four corners of the plot that is not used. In 
most situations the diagonal lattice provides better visuahzation than the 
natural lattice. This kind of non-diagonal lattice was used by Innui et al [15] 
with a slightly different shift operator. Our evolution equation, however, has 
the advantage of preserving the final probability distribution — except for the 
above mentioned rotation — ^with the same coin operator. 

In Fig. 4 we have results of two experiments of double slit with Hadamard 
walkers. The slits were placed exactly as in the previous example: one at 
(20, 6) and another at (20, —6). In the first simulation we carried out T = 100 
steps; in the second simulation, T — 800 steps. The simulation took about 



10 




Fig. 3. Probability distribution after a hundred steps of a Hadamard walker. Here, 
the shift operator is such that the mathematical lattice coincides with the physical 
lattice. 




20 40 60 80 100 120 140 160 180 200 200 400 600 800 1000 1200 1400 1600 

site number site number 



(a) Simulation with T = 100 steps (b) Simulation with T = 800 steps 
and screen along x = 60. and screen along x = 500. 

Fig. 4. Simulation of observation screens in the experiment of double slit. 

15min in the latter case. In the plots of Fig. 4 we have the pattern that would 
be observed in screens placed, respectively, along x — 60 and x — 500. In the 
latter case we note that the local minima are lower and the curve smoother. 



3.3 Simulation of detectors 



Detectors are described by a collection {Mm} of measurement operators which 
satisfy the completeness equation, J2rn ^'^LMm = I ■ The index m indicates the 
result of the measurement. Before state |\E') is measured the probability that 
result m occurs is given by p(m) = (^| M^^Mm \^), and after the measurement 
the state collapses to 1^') = }- — 1^). In qw2d we may define an arbitrary 

number of detectors specifying a list of coordinates (mi,ni), ■ • • , {mN,nj\[). 
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^ site number 

(a) Contour plot of the final proba- (b) Simulation of observation screen 
bility distribution 

Fig. 5. Results of a double-slit experiment with Grover walker. Both the wall and 
the observation screen are parallel to the secondary diagonal, and there is a detector 
near one of the slits. 

The measurement operators are 

N 

Mo = /4®/oo-E^i (11) 
1=1 

and 

Mj = /4 (g) \mi, rii) {rrii, n^l , (12) 

for 1 < i < A^. 

In Fig. 5 we have the plots of another double-slit experiment. In this simulation 
we use the Grover coin and place the wall parallel to the secondary diagonal. 
The wall goes from (—60, —100) to (100, 60). One slit goes from (13, —27) to 
(15, —25) and the other goes from (25, —15) to (27, —13). We have also placed 
a detector near the first slit, at (15, —27). The experiment was repeated ten 
times in order to take the average of the results. 

The positions of the detectors are defined in the main section of the input 
file, by using the DETECTORS keyword followed by the number of detectors and 
their respective coordinates. The number of repetitions of the experiment is 
defined by using the EXPERIMENTS keyword. Here we have an example on how 
to use these two keywords: 

DETECTORS 1 15 -27 
EXPERIMENTS 10 

The observation screen, also placed parallel to the secondary diagonal, goes 
from (20, —100) to (100, —20). The x axis in Fig. 5(b) is labeled sequentially, 
starting from the first point of the screen. We note in both plots that the 
interference pattern was asymmetric and also weaker on the detector's side. 
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3-4 Simulation of finite lattices 



The QWalk simulator can be used to simulate quantum walks on finite lattices. 
In this Section we study the example of a square lattice but QWalk also 
allows the definition of arbitrary boundaries. In the case addressed here, the 
boundaries were generated by breaking the links over a square the vertices of 
which are (-M,M), (M,M), (M, -M) and (-M, -M), as in [12]. Different 
values of M were investigated. The stationary distribution was approximated 
by running the simulation for 5000 steps. 

In order to calculate the total variation distance we must have the MIXTIME 
keyword in the main section of the input file, followed by the number of steps 
used to approximate the stationary distribution. This number should be higher 
than — or at least equal to — the number of steps that will be simulated, and 
the positions of permanent broken links must be correctly declared in order to 
define a closed region of the lattice. Since the number of steps simulated may 
be much greater than the size of the achievable lattice, we may improve the 
performance of the simulation by using the LATTSIZE keyword. This keyword 
declares that the lattice is only used from from x = —max to x = max and 
from y = —max toy = max, where max is the integer number passed as argu- 
ment of the LATTSIZE keyword. This option saves both memory and running 
time and should be used whenever the walk is restricted to a finite region of 
the lattice. The LATTSIZE keyword must be used after the STEPS keyword. 

For the previous example, taking M = 60, we may prepare an input file with 
the keywords 

MIXTIME 5000 
STEPS 2000 
LATTSIZE 59 

together with the keywords declaring the boundary, 

BEGINBL 

LINE -60 60 60 60 LINE 60 60 60 -60 
LINE 60 -60 -60 -60 LINE -60 -60 -60 60 
ENDBL 

When the MIXTIME option is used, the approximate stationary distribution TCapp 
is obtained in the beginning of the simulation. Afterwards the .sta output file 
keeps the total variation distance, at each step, from the average distribution 
to both the stationary and the unifornf^ distributions, so that the user may 
easily plot this information with an appropriate software such as gnuplot. 



^ In its present form, the simulator calculates the uniform distribution over all sites 
of the mathematical lattice, not taking into consideration the particular boundary. 
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(a) Total variation distance (b) Standard deviation 

Fig. 6. Total variation distance between the average distribution and the approx- 
imate stationary distribution as a function of time; and evolution of the standard 
deviation. Both plots refer to the Hadamard walk in the diagonal lattice and for 
different sizes of square boxes. 




Fig. 7. Stationary distribution approximated with 5000 steps and a box with 
M = 60. 

In Fig. 6(a) we have the total variation distance between the average distribu- 
tion and the approximate stationary distribution as a function of time, for the 
Hadamard walk in a finite lattice with rectangular boundaries, for different 
values of M. In Fig. 6(b) we have the evolution of the standard deviation of 
the Hadamard walk in the diagonal lattice for different values of M. We find 
in the .sta output file enough information to generate this kind of plot with 
just a few gnuplot commands, or with any other similar tool. We should note 
that these plots are also consistent with the results obtained by Oliveira et 
al [12]. 

In Fig. 7 we have the approximate stationary distribution for the Hadamard 
walk inside a box with M = 60, which is visually different from the uniform 
distribution. Moreover, one may also confirm with the output file from QWalk 
that the average distribution does not converge to the uniform distribution 
even for large times. 
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Fig. 8. Total variation distance as a function of time from the average distribution 
for the decoherent two-dimensional walk to both the uniform and the coherent 
stationary distributions. Two sources of noise are compared. The Hadamard coin 
was used and the stationary distribution was approximated with 7 • 10^ steps. 



3. 5 Simulation of decoherence 



QWalk allows the simulation of two-dimensional quantum walks with two dif- 
ferent sources of decoherence. The first of them is generated by measurements 
from randomly disposed detectors. It is also possible to simulate a unitary 
noise [16] which randomly opens links of the lattice. This noise model has 
been studied for the quantum walk on a line [9] and on a plane [8] . 

The measurement decoherence is better observed in the simulation of finite lat- 
tices. In order to declare the probability of measuring each site we use DTPROB 
keyword in the main section of the input file, followed by a numerical value. In 
Fig. 8 we compare two kinds of noises — measurements and broken links — with 
the same decoherence parameter p = 0.01 and the same box size M = 20. In 
both cases the experiment was repeated ten times with the Hadamard coin 
in order to take the average of the results. In the dotted lines we have the 
results for the measurement decoherence and in the full lines we have the re- 
sults for the unitary decoherence. In both cases we have the total variation 
distance, as a function of time, between the average distribution and both 
the uniform and the coherent stationary distributions. The stationary distri- 
bution was approximated with 7-10^ steps. We note that in both cases the 
average distribution initially approaches the coherent stationary distribution 
until ~ 1/p time steps, going to the uniform distribution after that. 

In Fig. 9 we have the result of a simulation of a Fourier walk with asymmetric 
unitary decoherence. The results are similar to those obtained by [8]. We 
declared a probability po = of broken links in the secondary diagonal, and 
a probability of pi = 0.2 of broken links in the main diagonal. This kind of 
simulation can be done with the BLPROB keyword, which must be followed by 
the probabilities of broken links in each direction — in the secondary and the 
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Fig. 9. Probability distribution after a hundred steps of a decoherent Fourier walker 
in the diagonal lattice. The probability of broken links was asymmetric, namely 
= in the secondary diagonal and pi = 0.2 in the main diagonal. 

in main diagonals, when the diagonal lattice is used; in the horizontal and in 
the vertical directions, when the natural lattice is used instead. 



3.6 Simulations in one- dimensional lattices 



The usage of qwld is analogous to that of qw2d, except for the absence of 
some keywords. Namely, qwld does not recognize SCREEN, BLPERMANENT and 
DETECTORS keywords and docs not recognize FOURIER and GROVER sub-options 
for coins and states. It docs not recognize DIAGONAL or NATURAL sub-options 
for LATTYPE either. On the other hand qwld works with three kinds of lat- 
tices, selected with LATTYPE keyword: the infinite line is selected by LINE 
sub-option; the finite one-dimensional lattice with reflecting boundaries is se- 
lected by SEGMENT sub-option; and the cycle is selected by CYCLE sub-option. 
Further information on the usage of qwld is provided together with the source 
(or binary) files. 

In Fig. 10 we have the probability distribution for a Hadamard walk in a one- 
dimensional infinite lattice [11,17,18,19]. We compare the coherent-evolution 
case with a decoherent one, where dccohcrence is introduced by random broken 
links. We could have introduced decoherence by means of random measure- 
ments as well. The probability of broken links in the example is p = 0.01. In 
both the coherent and in the decoherent cases we have executed T — 1000 
steps with qwld. In the latter case we have taken as result the average over 
a hundred independent experiments. We note in the figure that, due to the 
number of steps being much higher than 1/p, the classical behaviour of the 
quantum walker has already emerged [9] . 
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Fig. 10. Quantum walk on a one-dimensional lattice with broken links. 
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(a) Final probability distribution after (b) Stationary distribution approxi- 
T = 2 • 10^ steps mated with T = 10^ steps 

Fig. 11. Quantum walk on the cycle with a hundred sites. 

In Fig. 11 we have the result of a Hadamard walk on the cycle [10,11] with 
a hundred sites after T — 2 • 10^ time steps. Aharonov et al [13] proved 

that the quantum walk on a cycle with an odd number of sites mix to the 
uniform distribution. The same, however, does not hold in general for cycles 
with an even number of sites. This is clear in Fig. 11(b), where the stationary 
distribution was approximated with a large number of steps. One can also 
confirm this remark by using the information generated by QWalk to plot the 
total variation distance to both the uniform distribution and the approximate 
stationary distribution at each time step. The plot would show that the former 
does not converge to zero while the latter does. 



4 Conclusions 



In this work we initially reviewed the basics of the discrete quantum walks 
in infinite two-dimensional lattices with the possibility of broken links and 
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two different sliift operators. Tlie one-dimensional walk is a particular case of 
the one described. We presented then the QWalk simulator and described its 
usage. The QWalk is a quantum walk simulator for one- and two-dimensional 
lattices. We showed by means of examples how QWalk can be used to simulate 
double-slit experiments, walkers in finite lattices, detectors and two different 
kinds of decoherencc. Two kinds of lattices are available for two-dimensional 
walks and three kinds of lattices for one-dimensional walks. 

The simulations presented in this paper correspond to recent works in quantum 
walks, and have been performed with quite simple commands of the QWalk 

simulator. Apart from having simulated the available results of literature with 
great precision, the simulator is also an important tool for researchers to carry 
out new experiments. 

One of the most important potentialities of this simulator is the possibility 
of studying the walk with different boundaries, by defining permanent bro- 
ken links appropriately. It also makes possible to investigate the influence of 
detectors in the quantum walk, and the behaviour of the mixing time in dif- 
ferent situations. The simulator can be useful to simulate interesting physical 
settings, such as the double-slit experiment given as example. 

In future versions some improvements may be introduced in QWalk. The sin- 
tax of the input file, for instance, would be better if one could use numeri- 
cal constants instead of typing several digits when declaring custom coins or 
states. The treatment of the broken links may also be improved in order to 
simplify the definition of complex boundaries and to allow, for instance, the 
definition of time-dependent broken links, which would allow moving bound- 
ary conditions. It may also be useful for the researcher that the simulator saves 
output files with the wave-function in intermediate instants of the simulation. 
It should also be interesting to measure the amount of entanglement during 
the walk. 
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A Further options available in the simulator 



In this Appendix we describe the options that were not used in the examples. 
The CHECK keyword requests that some consistency checks be carried out dur- 
ing simulation. A sub-option is expected after this keyword: the STATEPROB 
sub-option declares that the unitarity of the state should be checked at each 
step of simulation; and the SYMMETRY sub-option, that the symmetry of the 
probability distribution is to be checked. Naturally, the latter sub-option 
should only be used when the symmetry of the probabihty distribution is 
suspected a priori, in order to confirm that suspicion. In two-dimensional 
simulations we may select two kinds of symmetry checking. The XSYMMETRY 
sub-option checks the symmetry around X clXlSj 1. e., checks if the probability 
at site {x,y) is the same as the probability at site (— a;, |/), for all x,y in the 
lattice. Analogously, we have the YSYMMETRY sub-options. 

In order to define a custom coin we must have a COIN CUSTOM command in 
the main section of the input file and must describe the coin on a separated 
section in the same file. This section must be delimited by BEGINCOIN and 
ENDCOIN keywords, and must contain all the entries of the coin matrix, from 
left to right and from top to bottom, with real and imaginary parts separated 
by a blank space. For instance, the code 



In two-dimensional simulations the description of the coin is analogous. 

In order to define a custom state we must have a STATE CUSTOM command 
in the main section of the input file and must describe the coin on a sepa- 
rated section in the same file. This section must be delimited by BEGINSTATE 
and ENDSTATE keywords, and must contain all the non-zero amplitudes of the 
initial state, with real and imaginary parts separated by a blank space. In 
one-dimensional simulations each of these amplitudes must be preceded by 
two integers: the first indicating the coin and the second indicating the posi- 
tion of the walker in that basis state. For instance, the code 



BEGINCOIN 
0.707106781186 0.0 
0.0 0.707106781186 

ENDCOIN 



0.0 0.707106781186 
0.707106781186 0.0 



defines, in one- dimensional walks, the coin 




(A.1) 
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BEGINSTATE 

1 0.0 1.0 
ENDSTATE 

defines, in one-dimensional walks, the initial state 

l*o)=^|l)|0). (A.2) 

In two-dimensional simulations the description of the state is analogous except 
that it requires two integers for the coin and two integers for the position of 
the walker. 

The SEED keyword manually sets a random seed for the random number gen- 
erator. The default is a number obtained from the system clock and usually 
the user should not change this. 

In two-dimensional walks the AFTERMEASURE keyword defines the number of 
iterations that will be carried out after a non-trivial measurement, i.e., after 
the result of a measurement corresponds to one of the detectors instead of its 
complement. 

LATTEXTRA keyword is used quite rarely. It defines an extra space to be re- 
served for the lattice in order to avoid access to invalid regions of memory 
during simulation. Its default value is 1 and normally this value should not be 
changed by the user. When this keyword is used together with LATTSIZE and 
STEPS there is a correct order that must be respected: first LATTEXTRA, then 
STEPS and finally LATTSIZE. Using LATTEXTRA keyword may cause QWalk to 
provide bad gnuplot scripts. In this case the user should fix the script, provid- 
ing appropriate ranges and "tics" . 
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